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Abstract 

We numerically study the metastable states of the 2d Potts model. Both of equilibrium and 
relaxation properties are investigated focusing on the finite size effect. The former is investigated 
by finding the free energy extremal point by the Wang-Landau sampling and the latter is done by 
observing the Metropolis dynamics after sudden heating. It is explicitly shown that with increasing 
system size the equilibrium spinodal temperature approaches the bistable temperature in a power- 
law and the size-dependence of the nucleation dynamics agrees with it. In addition, we perform 
finite size scaling of the free energy landscape at the bistable point. 
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I. INTRODUCTION 



Phase transition dynamics is ubiquitous in our everyday life. One of the most popular 
examples is boiling of water jl| . In spite of its familiarity, our quantitative understanding of 
such phenomena remains quite poor due to its strong nonequilibrium property. The most 
oopular theory of these first order transition dynamics is the classical nucleation theory 
2|, y| based on the droplet excitation picture, which fundamentally consider the cost of 
surface energy and gain of bulk energy. However it is beyond this theory how many nuclei 
are produced for unit volume in the macroscopic system 4^. It is necessary to consider 
multi-nucleation including the interaction between nuclei. 

Although nucleation is essentially nonequilibrium process, equilibrium property is useful 
as a starting point. First order transition, however, leaves more unclear points even in 
equilibrium in comparison with second order transition. One reason is that the former is very 
hard to investigate with conventional numerical simulations such as Monte-Cairo methods; 
in order to know fine metastable stricture we have to observe sufficiently large number of 
rare events, overcoming of free energy barrier. Extended ensemble methods, however, have 
made breakthrough to this problem 5ru|. In this article, we report the detailed analysis of 
equilibrium metastable states by the Wang-Landau sampling method and its relation to the 
transition dynamics. 



II. EQUILIBRIUM SPINODAL POINT 

To investigate the equilibrium met ast ability, we employ the Wang-Landau sampling 
method 0], which yields the density (number) of states as a function of macroscopic exten- 
sive quantities such as internal energy. This method is quite efficient to sample metastable 
states which appears with low probability in the Boltzmann distribution since this method 
realize a flat energy histogram. For 2d Potts model, we calculate density of states g(E) as 
a function of the internal energy, 

E= i 1 - <Wj) with ^ = 0, 1, • • • , g - 1. 

<i,j>Gn.n. 

The Helmholtz's free energy F((3) is written as 



where (3 = 1/ksT is inverse temperature and X denotes microscopic state. The extremal 
condition; dF((3; E)/dE = 0, gives the energy corresponding to the maximum or minimum 
probability in canonical ensemble at given j3. It is rewritten as 

P = ^fip- where S(E) = \ng(E). (1) 

When (3 is not monotonic function of E, there are multiple minima in canonical ensemble 
for given (3 indicating metastability and hysteresis. The things may be more clear to regard 
j3 in Eq. ([[]) as a response of thermodynamic function S(E) against the change of E in the 
microcanonical ensemble. In microcanonical ensemble, there are coexisting phase between 
two distinct pure phases. The phase transition between pure phase and coexisting phase 



has direct connection to t 
and lattice gas system 



re evaporation/condensation transition investigated in continuum 



Hill. 



The q-state Potts model on square lattice takes a first order transition between the 
paramagnetic and ferromagnetic phases for q > 4 at (3 C = ln(l + y/q). We investigate the 
system with q = 8 but believe that qualitative behavior does not depend on q. Periodic 
boundary condition is imposed for L x L square lattice. We perform parallel computation 
to treat large size system. The energy region is divided into some parts [l2] and each part 
is associated to individual thread. In order to enhance the relaxation and guarantee the 
ergodicity, we make overlap energy region for neighboring threads, where the exchange of 
spin configuration is allowed satisfying a detailed balance condition, 

W({X, Y} -» {Y, X}) _ 9l (E(Y)) 9J (E(X)) 
W({Y, X} -> {X, Y}) 9i (E(X)) 9j (E(Y)) ' 



in the same spirit with the parallel tempering method^]. Here, Qij Qj is a density of states 
calculated by z-th/j-th thread and Y} — > {Y, X}) means the transition probability 

from the compound state, X for the z-th thread and Y for the j-th thread, to its exchanged 
state. The density of states is calculated for the system with L = 32 — 2048 by using 128 
threads at most. 

The left panel of Fig. [1] shows relation between (3 and E with various system size. For 
given (3 around (3 C = 1.3424 • • • , there are two free energy minimal points and one maximum 
one between them. There are two spinodal points, (E^ pi , f3^ pi ) and (E^ pi , /3f pi ), at the end 
of low-energy ferromagnetic branch and high-energy paramagnetic branch as saddle node 
bifurcation points where energy minimal and maximal points meet to annihilate together. 
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FIG. 1: (left) inverse temperature as a function of internal energy calculated from difference, 
In g(E + 1) — In g(E), averaged over four individual simulations. The horizontal line indicates 
(3 = j3 c = ln(l + v8)- (right) size-dependence of the deviation of spinodal temperature from the 
first order transition temperature. 

The both spinodal points explicitly depend on the system size and approach f3 c with increas- 
ing L, i.e., the hysteresis region, /?/ pi < j3 < /3g pi , tends to disappear in the thermodynamic 
limit. The difference is plotted as a function of L in right panel of Fig. [TJ which indicates 



|/3 spi - p e \ oc L- 1 '" 



with 



i/ = i.4(i; 



(2) 



both for /3/ pi and /3f Di . 



spi 



In addition, the spinodal energy E f spi /N and E P S JN approaches E{/N « 0.40 and 



spi/ 



Ep/N w 0.89, respectively, with increasing L. These suggests a scaling behavior at the 
bistable point /3 C ; 



\E spi - £ c | oc L 



and furthermore 



F(p c ; E) = E- fc 1 lng(E) « F min + L**V ({E - E f c )/L dE ) , 

for El < E E v c . Hereafter we consider the scaling behavior around E[ but the same 
argument holds around E v c . Scaled barrier function V(x) rises for x C 1 and flat for x ^> 1. 
This leads that 

1 - L- {dE - dF) V'(E/L dE ) 



dE 




FIG. 2: finite size scaling of (3(E) (left) and F(f3 c ,E) (right). 

This reproduce Eq. fl2]) with <1e — cLf = \jv. We perform scaling as shown in Fig. [2] to 
estimate dp = 1-2(1) and dp = d>E — l/v = 0.5(1). Good collapsing of data is obtained for 
L > 256. Note that this finite size scaling does not hold for E < E£, where both internal 
energy and free energy are extensive, i.e., dp = dp = 1. The free energy well becomes more 
asymmetric as system size increasing. 

Naive energetics of droplet excitation suggests dE — d — 1 = 1. The slight enhancement 
of the excitation energy is presumably due to the fractal roughness of the interface. On the 
other hand, the dimension of free energy barrier dp = 0.5(1) is smaller than 1. This deviates 
from the behavior of macroscopic phase separation, dp — 1 [3] and there should be another 
crossover for larger scale, E — E{ ^> L 5 . 



III. NONEQUILIBRIUM RELAXATION, ESCAPE FROM METASTABLE 
STATE 

Next, we investigate transition dynamics from ferromagnetic state to para state at fixed 

nn 

temperature. The time evolution of the system is given by the Metropolis algorithm |3l. 1141] . 
The unit time is defined as the span in which all spins are updated one time in a sequential 
manner. All spins o~i are set to the same state named in the initial state at t — —1000. 
After the system is equibrated at (3 C , is suddenly lowered to aimed value at t — 0. 

We observe an order parameter (m) = (A r ~ 1 (^,-5 (Tii o) — q ,_1 )/(l — Q'^ 1 ), where (• • • ) means 
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FIG. 3: (left) time evolution of order parameter for L=256 and 1024. The data are averaged over 
16384/L samples, (right) temperature dependence of the time crossing m = 0.4 for various system 
size. The vertical lines denote the equilibrium spinodal temperatures for various system size. 

the average over independent random number realizations. Time evolution of (m) for L=256 
and 1024 is shown in the left panel of Fig. |3j For L=1024, there is almost no finite size 
effect in this data. The characteristic decay time tends to diverge as approaching C . We 
plot the time T0.4 at which (m(t)) crosses 0.40 as a function of (3 in right panel of Fig. |3j 
For sufficiently large L, size dependence is not observed and the divergence seems to be 
a little faster than a power-law type; r « (1 — (3/(3 c )~ z,y . This implies that the dynamics 
is of thermal activation type, typically r oc exp [const. /(/3 C — /?)], although it is difficult to 
speculate its analytic expression. 

On the other hand, finite size effect appears approaching (3 C the earlier, the smaller L 
becomes. The inverse temperature above which T0.4 deviates from large size behavior agrees 
well with the finite-size spinodal point Pg pi (L) obtained by equilibrium simulations. Initially 
T0.4 overcomes the large size behavior, which presumably due to the backward transition 
from para phase to ferro phase, while T0.4 becomes smaller than large size behavior and 
saturates to a certain value increasing with L. 
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IV. SUMMARY 



In conclusion, we calculate the size-dependent spinodal points of the 8-state Potts model 
and perform finite size scaling to obtain fractal exponent for free energy barrier dp ~ 0.5 
and internal energy dg ~ 1.2. These exponents can be compared with (d 2 — d)/(d + l) = 2/3 
and d 2 /(d+ 1) =4/3 in Ref. 10], where evaporation/condensation transition is discussed 
based on the Ising model in micromagnetic ensemble, by relating /3 spi — j3 c and (E — E c ) /N 

~ (2) 

to magnetic field H\ and magnetization m t — m coex , respectively. 

Additionally it is confirmed that size dependent spinodal temperature also gives the 
dynamical criterion in transition dynamics, which suggests some connection between equi- 
librium metastability and the nonequilibrium relaxation. Equation [2] implies that there is a 
diverging length scale in bulk 

* I/ 3 - 

with v PS 1.4. It is important to note that such critical like behavior is observed in the 
energy region of coexisting phase, which is not realized in the canonical ensemble. It will 
be very interesting if we can find any relation between these unphysical equilibrium state 
and the nonequilibrium transient state with corresponding energy. The exponent v ps 1.4 
is larger than 1 which is expected for critical nucleus radius and may mean the average 
distance between critical nuclei. The advanced study is required for the thermally activated 
nuclear growth, where typical distance of nuclei obtained here will give a useful hint. 
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